Reversal of the circulation of a vortex by quantum tunneling in trapped Bose systems 
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We study the quantum dynamics of a model for a vortex in a Bose gas with repulsive interactions 
in an anisotropic, harmonic trap. By solving the Schrodinger equation numerically, we show that the 
circulation of the vortex can undergo periodic reversals by quantum-mechanical tunneling. With 
increasing interaction strength or particle number, vortices become increasingly stable, and the 
period for reversals increases. Tunneling between vortex and antivortex states is shown to be 
described to a good approximation by a superposition of vortex and antivortex states (a Schrodinger 
cat state), rather than the mean- field state, and we derive an analytical expression for the oscillation 
period. The problem is shown to be equivalent to that of the two-site Bose Hubbard model with 
attractive interactions. 
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Over the past decade, atomic Bosc-Einstein conden- 
sates have provided unprecedented opportunities for 
studying in detail the properties of quantized vortices, 
both theoretically [l[ and experimentally [2j . Most theo- 
retical work has been based on the use of the mean-field, 
Gross-Pitaevskii (GP) approximation Q, and an inter- 
esting prediction by Garcfa-Ripoll et al. Q within this 
approach is that a rotating Bose-Einstein condensate in a 
non-rotating anisotropic harmonic trap can undergo pe- 
riodic reversals of the sign of the vorticity if the initial 
energy of the system is sufficient to overcome the energy 
barrier between vortex states with opposite circulation. 

Recent experimental developments make it possible to 
realize few body systems trapped on an optical lattice 
[H, Q . The behavior of such small systems can be quite 
different from what is predicted from a classical treat- 
ment based on the GP approximation. In this paper we 
consider the problem of stability of a vortex quantum- 
mechanically, and show by solving the Schrodinger equa- 
tion numerically that for energies below the barrier, re- 
versals of the vorticity can occur by tunneling. We find 
that the wave function during reversals resembles more 
closely a quantum superposition of states (a Schrodinger 
cat state) than a mean-field one, and we derive an analyt- 
ical expression for the rate of reversals that agrees well 
with the numerical data. Mathematically, the problem 
is equivalent to that of particles with attractive interac- 
tions in a double-well potential, which has been studied 
previously @, HI- 

Consider N identical bosons of mass m in an 
anisotropic, harmonic two-dimensional trap, with trap 
frequencies denoted by u x and u> y , and for definiteness we 
shall assume that lu x > oj v . Following Ref . [2] we shall as- 
sume that the only oscillator levels occupied are the first 
excited states of the oscillators @, corresponding to the 
wave functions ip x = C (x/d x ) exp (— x 2 /2d x — y 2 /2d 2 ), 

where C — \/2/ (nd x d y ) and d x = \Fhpn 



expression for ip y . The many-body Hamiltonian is 
H = Ho + Hi n t 

= E 
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where e, = hw+hoji (i = x, y) andw = (u) x +u) v ) /2, while 
cj creates and c, destroys a particle in the state tpi [l(|. 
We shall take the interaction to have the contact form 
(r, r'|V|r, r') = 52D^ 2 (r — r'), where g 2 D is an effective 
two-dimensional interaction strength, which we shall take 
to be positive [111 ]. 

For orientation we first describe the GP approach in 
which all particles are assumed to be in the same single- 
particle state. Thus the many-body state may be written 



|A0, An) ^^L= {n^ci+nl^cle-^Y |0>. (2) 
The quantities n x and n y are the occupation probabilities 



of the two states, and 



= 1 and Art = 



The 



phase difference between the two components is denoted 



by A0 



where <p x and 4> y are the phases of the 



two states, and |0) is the vacuum state. The GP energy 
functional is [12] 

|=/ d 2 r||l|V,M 2 +^(^^+^ 2 )|^| 2 + ^|^|^ 



(3) 

where ip is the condensate wave function, ip = a x tp x + 
a y t(j y = loa-le^Va: + \a y \e^^ y with la^ 2 + |a y | 2 = N. 
Evaluation of Eq. ([3]) for this wave function leads to 



£ = 



E 1 Aw . 7 r „ ,, .,, . ». , 



(4) 



where 7 = ^{)^Ju) x ui y /w 2 , with 70 = Na s /Z [111 ], and 
Alo = lu x — Lo y . In Eq. (j4|) we have omitted the con- 
and a similar tribution independent of An and Atfi [l3j . One sees 
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FIG. 1: (Color online) Landscape of the energy per particle £ 
for T = (a), 0.8 (b), 4.0 (c), and 8.0 (d). Purple regions (the 
darker area for An < 0) correspond to lower £ and red ones 
(the darker area for An ~ 1) to higher £. For a given panel, 
the contour lines are equally spaced in £ , but the spacing 
varies from panel to panel. 



that there is one dimensionless parameter in the prob- 
lem, T = 7u>/ Aui [3. 

In Fig.[Tl we show contours of £ as a function of A(f> and 
An for several values of the interaction strength. States 
with A<j) = — 7r/2 are vortex- like with positive circulation, 
while those with A(f> = n/2 correspond to an antivortex, 
a vortex with negative circulation. States with A<p = 0, 7r 
have nodal lines which, because of the larger mean square 
density, give rise to a larger repulsive energy. Allowed 
motions correspond to contours of constant energy. For 
small T [Figs.[lja) and (b)], all contours are open, and the 
allowed motions correspond to oscillations between vor- 
tex states and antivortex ones. For T > 1, global minima 
of the energy develop on the lines A<j) = —ir/2 and tt/2, 
and there are closed contours surrounding the minima. 
These correspond to motions in which a vortex line os- 
cillates without reversals of the circulation. With further 
increase of T, the closed contours occupy an increasing 
fraction of the area [see Figs. QJc) and (d)] and an en- 
ergy barrier with height ~ jhu) ~ Tiu)Na s /Z per particle 
grows between the vortex and antivortex states, i.e., large 
r stabilizes the vortex and antivortex states. Classically, 
for a system having energy less than the minimum value 
on the line Acj) = 0, the circulation of a vortex cannot 
change sign. 

However, in quantum mechanics the circulation can re- 
verse by tunneling. The interaction Hamiltonian is given 
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FIG. 2: (Color online) Time evolution of I for P = 4.0 and 
TV = 2 (a), 3 (b), and 5 (c). Panel (d) is the correspond- 
ing energy landscape showing that the vortex and antivortex 
states lie on closed orbits, and the arrow indicates a possible 
trajectory for tunneling between these states. 



by #int = {jhuj/4:N)h int with 

^int = 



V V 

+3(clc' x c x c x + r t 
= 3N 2 -2N-1 2 



Ac x CyCyC x 



c y c yCyCy) 



(5) 



Here N 

1 = 



c t c x + 



|0), corresponding to | — 7r/2,0) for the par- 



0.01. Preliminary results were 



is the total number operator and 

The expectation value L of the z- 
component of the angular momentum operator in a state 
| "J) containing only the two single-particle states ip x and 
ijjy is given by L = Ti^J Co 2 / u> x Lo y {^\l\^!) = Ti^Jui 2 /uj x uiy I. 
The system has N + 1 Fock states \N x ,N y ) = \N X ,N- 
N x ) with N x — 0, 1, . . . , N, where N x and N y are the 
numbers of particles occupying the single particle states 
tj) x and ij)y, respectively. Using the above many-body 
Hamiltonian, we have followed the time evolution of the 
system, starting from the vortex state proportional to 

ticular choice AujJuj 
reported in Ref. [15J. 

Classically, reversals of the circulation with this initial 
state are impossible if T > 2, which ensures that the 
state lies on a closed orbit. However, we see, e.g., in Fig. 
[5]for r = 4.0, that reversals of the angular momentum do 
occur. The angular momentum oscillates with a period 
T ~ 2.41 T for N = 2, T ~ 5.73 T for N = 3, and T ~ 
37.7 To for N — 5, where To = 2tt/Auj is the period in the 
absence of interactions. The rapid wiggles in Figs. 15(a)- 
(c) are due to oscillations of internal degrees of freedom of 
the vortex, which correspond to motion on closed energy 
contours in the mean-field approach. 

In Fig. [3(a), we plot the oscillation period T as a func- 
tion of N for r = 4.0 and 8.0. The results do not change 
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FIG. 3: (Color online) Oscillation period T as a function of 
N at fixed interaction strength T (a) and as a function of V at 
fixed N (b). The thick solid lines for V > 2 are calculated from 
the analytic formula ((8j). The lines connecting the data points 
in (b) are to guide the eye. The period T for the data points 
at r = 1.6 in (b) is hard to define and has an uncertainty of 
order 0.1 T . 



qualitatively for a more general model in which six single- 
particle states, the ground state and the three lowest d- 
like states in addition to the two p-like states, are taken 
into account [l(| • We see clearly that T increases almost 
exponentially with N. This tendency is consistent with 
the fact that, for these values of T, the vortex and an- 
tivortex states are stable classically. We also see that T 
increases with T. This is because for T ^> 1 the Hamil- 
tonian is dominated by H lrA and therefore level spac- 
ings are proportional to 7. Thus mixing of states by the 
anisotropy becomes less important for large T. For Y < 1, 
the oscillation period is approximately To, as one would 
expect from the fact that under these conditions rever- 
sals can occur in classical mean-field theory for all initial 
conditions. For T -C 1, we see repetitive collapses and 
revivals of L with a period ~ 2ttN/ju> [13] • 

In Fig. [3jb), we plot T as a function of T for fixed 
N and one sees the different behaviors for T < 1 and 
r > 1. The oscillations in the former region are well 
described by mean-field theory and T is almost constant 
with a value ~ Tq. For large T, vortex-antivortex oscil- 
lation occurs for the given initial conditions only in the 
quantum-mechanical calculation, and T shows a power- 
law dependence on T. 

To understand the behavior of the tunneling time, 
we have investigated the state |^(t))- Figure [4] shows 
its overlap K^IA^opt, 0)| with the mean-field state ((2]), 
where A0 O pt is chosen to maximize the overlap. The ex- 
pectation value of I, which is proportional to that of the 
angular momentum, is also plotted, and it is denoted by 
[for |¥(t)) and f opt for |A0 opt ,O). For T < 1, |*(t)) is 
well described by the mean-field state, and A<j> op t changes 
continuously. For Tq = 2.4 and 8.0, A</) opt jumps essen- 
tially discontinuously between approximately —tt/2 and 
7r/2, corresponding to the vortex and antivortex states. 
This indicates that the main components of the state cor- 
respond to either all particles being in the vortex state 
or all of them in the antivortex state. A much better 
approximation for the wave function for r > 1 is the 
Schrodinger cat state consisting of a superposition of the 
states in which all particles are in the vortex state or all 



I Cat; 9) = cos - 
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7T 
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(6) 



In Fig. 03 we plot the overlap K^ICat; O pt)| for the same 
|\&) as in Fig. 0Jd) .18j. Similarly, 8 opt is determined by 
maximizing |(\&|Cat; O pt)|- Here L opt is the expectation 
value of the angular momentum in the state |Cat; opt ). It 
is remarkable that 8 opt changes almost continuously and 
the overlap is close to unity (its mean value is ~ 0.97 for 
this case). 
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FIG. 4: (Color online) Overlap |(*|A^ op t, 0)| of the wave 
function with the optimized phase state for N — 5 and T = 
(a), 1.6 (b), 2.4 (c), and 8.0 (d). 
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FIG. 5: (Color online) Overlap | (^ICat; O pt)| of the numerical 
solution of the wave function with the cat state for N = 5 and 
T = 8.0. 

The behavior of the wave function for T ^> 1 
is most conveniently analysed in terms of 
states in which the interaction energy is di- 
agonal. These may be written as |Z) = 

{[(N + 0/2]! [(N - l)/2]\y 1/2 (c^)( 7V+/ )/ 2 (c t _)( Ar - / )/ 2 |0) 
where the operators = (cj. + icL)/y/2 and 

c_ = (cj. — ic^ y )/\/2 create particles in a vortex or 
an antivortex state, respectively. The operator / in the 
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subspace of states we are considering is c i + c + 
and the Hamiltonian reduces to 



H=—-(clc-+clc+) 



~2N 



N 



(c + C + C+C+ + C_C_C-C-) 

(7) 



This is equivalent to that of the Bose Hubbard model 
for two sites if + and — are regarded as site labels, 
with hopping matrix element —TiAuj/2 and on-site 
interaction —jhui/N @, Q. The anisotropy of the trap 
couples states of different / according to the Hamiltonian 
W = (TlAuj/2)(c+c- + c l _c + ), which we may treat 
perturbatively in this regime. The wave function is 
dominated by the vortex and antivortex states, and 
components from other states are suppressed by the 
large interaction energy for these states. Because of 
the anisotropy, the vortex and antivortex states are 
not energy eigenstates and the leading contributions 
to the mixing of the states for large T may be cal- 
culated by perturbation theory. Since the anisotropy 
couples only states in which I differs by 2, it is of Nth 
order in Aw. The leading contribution to the matrix el- 
ement mixing the vortex and antivortex states is AE/2 = 
H'_ N _ N+2 (AE^ N+2 y 1 ■ ■ ■ H' N ^ 4N _ 2 (AE N ^ 2 )~ 1 H' N _ 2 
where H' u+2 = {l\H'\l + 2) and AEi = E t -E N , with 
Ei = -("/hij/4N)l 2 . Here AE is the splitting of the two 
lowest states. For T 1 one finds [?J 



AE = MAluIN 



a(N) 



N-l 



(8) 



Here a(N) = TV[(TV - l)!]- 1 /(A f -i)/2 is a function that 
depends weakly on TV: a(2) — 1 and a(oo) — e/2. The 
oscillation period is given by T = 2irfr/ 'AE. In the strong 
interaction regime, AE given by Eq. ([5]) accounts very 
well for the numerical results, as is shown in Fig. [3] We 
remark that rotation of the trap will suppress the vortex- 
antivortex oscillations if 2Nhil > AE, where f2 is the 
rotational angular velocity. 



In conclusion, we have calculated rates for tunnel- 
ing between vortex and antivortex states in a simplified 
model without making a mean-field approximation and 
have obtained an analytical result in the limit of small 
tunneling. We have demonstrated that, in this limit, the 
wave function is much better described as a Schrodinger 
cat state that is a superposition of a state in which all 
particles are in the vortex state and one in which all par- 
ticles are in the antivortex state. We have also shown 
that the problem is equivalent to a two-site Bose Hub- 
bard model with attractive interparticle interactions. On 
the experimental side, recent advances in creating few- 
body systems trapped on an optical lattice 0, [g] offer 
the possibility of observing experimentally the effects we 
predict. 
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